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CN ' Abstract 

o. 

The energy evolution of average multiphcities of quark and gluon jets 
^^ I is studied in perturbative QCD. Higher order (3NL0) terms in the 

perturbative expansion of equations for the generating functions are 






found. First and second derivatives of average muhiphcities are calcu- 
lated. The mean multiplicity of gluon jets is larger than that of quark 
jets and evolves more rapidly with energy. It is shown which quantities 
^S^ ■ are most sensitive to higher order perturbative and nonperturbative 

corrections. We define the energy regions where the corrections to 
different quantities are important. The latest experimental data are 
discussed. 
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1 Introduction 

The progress in experimental studies of separated quark and gluon jets is 
very impressive [|l|. Their separation has become possible due to new meth- 
ods and high statistics at the Z^ resonance, which also allows the selection 
of jets with different sub-energies. Therefore the detailed study of the energy 
evolution of such properties of jets as their average multiplicities and widths 
has also become possible. It is well known that the average multiplicities of 
quark and gluon jets increase quite rapidly with energy but their ratio has a 
much slower dependence and tends to a constant in asymptotics. The energy 



dependence of this ratio and of the ratio of the derivatives of mean multi- 
phcities is of much debate nowadays, with new experimental data appearing 
rather regularly. The experimental procedure of the jet separation, the scale 
choice, and the manner in which soft particles are ascribed to jets, is often 
disputable, however. Therefore, in the following, we do not insist upon a 
direct quantitative comparison of data with theory but rather describe the 
correspondence between our present experimental and theoretical findings 
with the purpose of finding quantities sensitive to areas of disagreement. 

On the theoretical side, perturbative QCD (pQCD) provides quite definite 
predictions which can be confronted to experiment. The asymptotical value 
of the ratio of average multiplicities in gluon and quark jets equals 2.25 [@], 
much larger than its experimental values which range from about 1.04 at 
the comparatively low energies of the T resonance to about 1.5 at the Z° 
resonance |^, Q. The next-to- leading order corrections to this ratio r (we call 
them NLOr) reduce it from its asymptotical value by about 10% at the Z° 
energy |^, 0, |H- The next order NNLOr terms diminish it further, although 



calculations by different methods |^ |T0|] lead to different results as we discuss 



later. We have calculated the next-to-next-to-next-to-leading order (3NL0r) 
terms using QCD equations for the generating functions. The 3NL0r terms 
are small but nonetheless provide an improved description of data, as shown 
below. 

The computer solution of QCD equations for the generating functions has 
shown even better agreement with experiment for this ratio [|lll , and excellent 



agreement for higher moments of the multiplicity distributions as well [12 
Being perfect at the Z^ energy, the agreement for the ratio is less good at the 
T resonance where the theoretical curve is 15-20% above the experimental 
result (see Fig. 2 in Jill). -'■^ other words, the theoretically predicted slope of 



the ratio of multiplicities in gluon and quark jets is noticeably smaller than 
its experimental value. Nonetheless, there has been a steady convergence of 
theory and experiment as more accurate calculations and more sophisticated 
measurements have become available. Moreover, it is even surprising that 
any agreement is achieved given that the expansion parameter is extremely 
large (about 0.5) at present energies. 

Recently, experimental data about the ratio of slopes of average multiplic- 
ities in gluon and quark jets have been reported |T3| and the corresponding 



analytical expressions have been calculated ^^. The importance of studying 



the slopes stems from the fact that some of them are extremely sensitive to 



higher order perturbative corrections, and to nonperturbative terms at avail- 
able energies, as is discussed below. Thus they provide an opportunity to 
learn more about the structure of the perturbative series from experiment. 
The experimental group |T^ claims that it is possible to simultaneously 



describe the scale dependence of multiplicities in gluon and quark jets only if 
nonperturbative effects contribute constant terms to the multiplicities. We 
examine this issue. The success of the computer solution suggests an al- 
ternative possibility, namely that a combination of nonperturbative terms 
which rapidly decrease with increasing energy and higher order perturbative 
terms can describe the data, effectively replacing the constant offsets assumed 
in [|l3l . To verify this, one needs analytical calculations of the corrections for 



experimental quantities most sensitive to them. With this aim in mind, we 
consider analytical expressions for the above ratios, as well as for the ratio 
of the second derivatives (curvatures) of average multiplicities. These last 
ratios have not yet been measured. The general problem of perturbative and 
nonperturbative terms in the QCD equations for the generating functions is 
also addressed using simple analytical estimates. These estimates show that, 
at asymptotically high energies, the role of the nonperturbative region is not 
as important as the higher order perturbative terms but that it becomes more 
important at lower energies. The location of the borderline depends on the 
quantity under consideration. 

For the comparison of theory with data, we utilize the principle of lo- 
cal parton hadron duality (LPHD) [l^, i.e. the hypothesis that a hadron 
distribution (experiment) can be related to a parton distribution (theory) 
by a simple normalization constant. LPHD has been shown to be reliable 
for distributions dominated by soft particles, such as the multiplicity related 



quantities studied here (for a fuller discussion, see [0). 



2 Evolution of mean multiplicities in pQCD 

In perturbative QCD, the general approach to studying multiplicity distribu- 
tions is formulated in the framework of equations for generating functions (for 
reviews see 



|17| , p!q|). From the generating functions, one derives equations 
for average multiplicities and, in general, for any moment of the multiplic- 
ity distributions |jlO|, [11, ^H]. In particular, the equations for the average 



multiplicities of gluon and quark jets are 

{naiy))' = f dx-f^[K^{x){{nG{y + Inx)) + {ndy + ln(l - x))) - {nciy))) 



+nfK^{x){{np{y + \nx)) + {np{y + ln(l - x))) - {nciy)))], (1) 
{npiy))' = I dx-ilKf{x){{nG{y + hix)) + {nF{y + Hl-x)))-{np{y))). (2) 



From these equations one can predict the energy evolution of the ratio of mul- 
tiplicities r and the QCD anomalous dimension 7 (the slope of the logarithm 
of the average multiplicity in a gluon jet), defined as 

r = 7 — r , 7 = -, — r = In n^ . (3) 

{nF) {nc) 

Here, the prime denotes the derivative with respect to the evolution param- 
eter y = ln(p6/(5o)) P and are the momentum and initial angular spread 
of the jet related to the parton virtuality Q = pQ/2 (for small angles), 
(5o=const, K^s are the well known splitting functions, (no) and (np) are 
the average multiplicities in gluon and quark jets, {no) is the slope of (nc), 
and Hf is the number of active flavors. The perturbative expansions of 7 and 
r are 

7 = 7o(l - ai7o - 0270 " «37o) + 0(7o), (4) 

r = ro(l - ri7o - r27o - ^37o) + 0(7o), (5) 



where 70 = j2Ncas/'n, as is the strong coupling constant 

27r 



Poy 






+ 0{y-% (6) 



/3o = (llA^c - 2n/)/3, pi = [17N^ - nf{5Nc + 3Ci.)]/3,ro = NjCp, and in 
QCD A^c = 3 is the number of colors and Cp = 4/3. 

The limits of integration in eqs. (0), (§) are generally chosen either to be 
and 1 or e^^ and 1 — e^^. This difference, being negligible at high energies 
y, becomes more important at low energies. Moreover, it is of physical sig- 
nificance. With limits of e~^ and 1 — e~^, the partonic cascade terminates at 
virtuality Qo as seen from the arguments of the multiplicities in the integrals. 
With limits of and 1, the cascade extends into the nonperturbative region 



with virtualities smaller than Qq/2, i.e. Qi ~ xpQ/2 and (^2 ~ (1 — x)pQ/2, 
where the very notion of partons becomes scarcely applicable. This region 
contributes terms of the order of e~^, power-suppressed in energy. It is not 
clear whether the equations and the LPHD hypothesis are valid down to the 
cutoff Qq only or if the nonperturbative region can be included as well. 

Nonetheless, the purely perturbative expansions (^) and (H) with constant 
coefficients aj,rj and energy-dependent 70 are apphcable only in the case of 
limits and 1. This is quite natural for infrared-safe quantities like mean 
multiplicities. According to (|), each derivative gives rise to the factor 7. 
Therefore, a Taylor series expansion of eqs. (|1|) and (0) is equivalent [|19 



to a perturbative series. Terms of the same order should be equal on both 
sides. From this, one obtains the coefficients ai,ri. The analytical formulas 
for Oj, Ti with i < 3 are given in Appendix 1. For i <2, they were presented 
earlier |^, [7], ^. We present them here using a more general notation. The 
numerical values of aj,rj for different numbers of active flavors Uf and in 
SUSY QCD P are given in Table |l|. 
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0.185 


0.426 


0.189 


0.280 


- 0.379 


0.209 
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0.191 


0.468 


0.080 


0.297 


- 0.339 


0.162 


5 


0.198 


0.510 


-0.041 


0.314 


- 0.301 


0.112 


S 











0.188 


-0.190 


-0.130 



Table 1: Numerical values of the perturbative corrections up to order 7q 
for the multiplicity ratio r and the QCD anomalous dimension 7, based on 
integration limits of the generating functions from to 1 (see text); Uf is the 
number of active quark flavors while S refers to SUSY QCD. 



At the Z^ resonance, for which the energy of an unbiased gluon jet is 
about 40 GeV |^, Q , the subsequent terms diminish the value of r compared 
with its asymptotic value rg = 2.25 by approximately 10%, 13%, and 1% for 
Hf = 4. The smallness of r^ is a good indication of the convergence of the 
series at this energy. 

Note that we have used the traditional method (see |T^) of considering r 
and 7 to be the main ingredients of the calculation and that we thus employ 



the conventional definition for tlie order of the corrections. This method 
introduces some asymmetry between quark and gluon jets. A more natural 
way to structure the calculation might be to define the anomalous dimensions 
of gluon and quark jets symmetrically, as 

7 = ^ = Nr^c)r; 7^ = ^ = [lnK)]' (7) 

and to call the subsequent corrections next-to- leading-order (NLO), etc. 
With the conventional approach used here, it is necessary to differentiate 
between corrections to r (which we call NLOr, etc.) and those to 7. For 
example, the NLOr-correction given by ri contributes to ''jp only at NNLO 
and higher orders (for more details see the next section). The differences 
between the traditional and symmetric methods are not important, however, 
as long as the corrections themselves are small. 

The behavior of the ratio r in various approximations (DLA=leading 
order, NLOr, NNLOr, 3NL0r) is shown in Fig. 1 for nj = 4. The abscissa 
variable is k = pQ ~ 2i?sin(9/2), where y ~ hi^n/Qo). Here Qo is chosen to 
depend on the number of active flavors in the proportions used previously pO 



in order that as agree with experiment. A comparison of these results to 
data is given in section ^. 

Now we would like to return to the issue of the limits of integration 
mentioned above. One can analytically estimate how the coefficients a^, r^ are 
modified if the range of integration becomes somewhat narrower, e.g. from 
e to 1 — £ (e =const ^1). In that case one easily calculates the corrections 
to aj,rj. These are given in the integrals Vj presented in Appendix 1. In 
particular, using the formulas of Appendix 1 for e = 0.1, nj = 4, one 
obtains rather large contributions from the nonperturbative region: ai^ = 
— 0.36ai, r^^ = — 0.56ri. However e~^ = 0.1 only at the T-resonance and 
is much smaller at the Z^. Therefore, narrowing the range of integration 
reduces the energy dependence. The 3NL0 result for r based on the narrower 
range of integration is shown by the curve labeled r{e) in Fig. 1. Surely, it 
is not entirely satisfactory to replace e~^ in the limits of integration by a 
constant e because if e depends on y the coefficients ai,ri become energy 
dependent and their derivatives need to be considered as well. Nonetheless 
such a procedure can be used to obtain an order-of-magnitude estimate. 
To obtain better accuracy in the energy region from the T to the Z^, the 
computer solution of eqs. (1), (2) with modified (in the nonperturbative 
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region) coupling strength, higher order terms in the kernels and more strict 
conservation laws (in particular, of transverse momenta) is necessary. 

3 Slopes as a stumbling-block of pQCD 

We have shown how the perturbative series can be used to evaluate 7 and 
r, and where nonperturbative corrections become important. In this section, 
we show that the perturbative approach should work well for the ratio of 
the slopes of the multiplicities but that it is much less reliable to use low 
order perturbative estimates, even at the Z°-energy, for such quantities as 
the slope of r or the ratio of slopes of logarithms of multiplicities (we shall 
call them the logarithmic slopes). Much larger energies are needed for that. 
Thus the values of r and/or of the logarithmic slopes can be used to verify 
the structure of the perturbative expansion. 

The slope r' is extremely sensitive to higher order perturbative correc- 
tions. The role of higher order corrections is increased here compared to r 
because each nth order term proportional to 7q gets an additional factor n 
in front of it when differentiated, the main constant term disappears, and 
the large value of the ratio r2/ri becomes crucial: 



5rori7o 



2r27o , /3r3 „ \ o , ^/ sn 



l + ^+(^^ + 5i)7o^ + 0(7o^) 



(8) 



where the relation 7q ^ —B'~^l{\ + -Bi7o) ^^^ been used with B = Po/8Nc = 
1/2C^; Bi = (3i/4:Nc(3o. The factor in front of the bracket is very small even 
at present energies: BtqTi k, 0.156 and 70 ~ 0.5. Nonetheless, the numerical 
estimate of r is unreliable because of the expression in the brackets. Each 
differentiation leads to a factor as (or 7q), i.e. to terms of higher order. For 
the values of ri, r2, r^ given in Table 1 (rij = 4), one estimates ^r^jrx ^ 4.9, 
(Srs/ri) + i?i ^ 1.5 (see eq. ( |42| ) and note the aforementioned factors of 
2 and 3 in front of r2 and rs, correspondingly). The simplest correction, 
proportional to 70, is more than twice as large as 1 at present energies, while 
the next correction, proportional to 7q, is about 0.4. Therefore even higher 
order terms should be calculated before the perturbative result for r' can be 
considered to be reliable. 

The slope r' equals for a fixed coupling constant. In the case of the 
running coupling constant, it evolves rapidly with y in the lowest perturbative 
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approximations (|^) (this is also seen from the Figs, of |Tl[), and ranges 
according to (|D from 0.06 at the Z° to 0.25 at the T. 

Let us consider the ratios of slopes r*^^^ and curvatures A"^^ defined as 

Their ratios to r can be written as 

Pi = -?!) = 1 - -, (10) 

r Z'yrr + rr — 2r , 

^2 = ^ = 1 (^2^y)^2 • (11) 

Since r oc 7q, the asymptotical values of r, r^^\ r*^^^ coincide and should 
equal 2.25. Moreover, the values of r, r^^\ r^"^^ coincide in NLOr. This im- 
plies that the mean multiplicities in gluon and quark jets and their derivatives 
differ in NLO by a constant factor tq only. First preasymptotical corrections 
are very small. They are of the order of 7q with a small factor in front and 
contribute about 2-4% at present energies: 

pi = 1 - 5^7^ ^ 1 - 0.077o', (12) 

P2 = 1 - 25ri72 ^ 1 - 0.147^. (13) 

However, this ideal perturbative situation does not persist when the next 
terms are calculated (see Appendix 2). For example, for pi their numerical 
values are 

pi = 1 - 0.077o'(l + 5.3870 + 4.1472) (n/ = 4). (14) 

All the correction terms are much smaller than 1. However, the numerical 
factors inside the brackets are so large that at present values of 70 ~ 0.5 the 
subsequent terms are larger than the first one and the sum of the series is 
unknown. One can sum the series using the simplest Pade-approximant (i.e. 
by assuming the steady decrease of the unaccounted terms in proportions 
determined by the ratio of two last calculated terms), resulting in 



pi = 1 - 0.077o' 



5.3870 
- 0.7870 



(15) 



In contrast to pi, the perturbative corrections to the ratio of slopes r^^^ 
are small. The lowest order 0(70) correction to r'-^-' is the same as for r but 
higher order corrections are smaller because they are negative both in r and 
Pi which define r^^^ = r/pi. This explains why experimental values of r*^^'' are 
similar to values of r calculated in the NLOr-approximation (see section ^), 
whereas experimental values of r are still about 25% lower than the NLOr- 
prediction. Yet more substantial cancelations of higher order terms should 
occur for r^^\ It would be interesting to check this prediction experimentally. 

The value of pi also determines the ratio of the anomalous dimensions, 
i.e. of the slope of the logarithm of the average multiplicity of a quark jet 
{■jf) to that of a gluon jet (7): 

r 

lF = Pil = 1 • (16) 

r 

The logarithmic slopes of quark and gluon jets are equal in NLO. They differ 
in higher orders in such a manner that 7f < 7 since both r and r are 
positive. Our failure to obtain a precise estimate of pi implies that we can 
not perturbatively evaluate '-)p either. The ratio of logarithmic slopes pi is 
more sensitive to higher order perturbative terms than the ratio of slopes 
r'^^\ similar to r . 

An interesting feature of eq. ([T^) is that pi contains terms up to Oipf^) 
whereas the r in the numerator of (|1^) is known to 0(7o) only. This situation 
is related to the fact that r'/7 ~ 0(7q) and allows the perturbative expression 
for '-)p to be written as: 

7^ = 7o [ 1 - ai7o - (^2 + Bn)-il - («3 + 2Sr2 + BrD^^l 

-(a4 + S(3r3 + 3r2ri + B^r^ + rl))^^ ] (17) 

if the higher order 4NL0 term —0479 i^ added to 7 inside the brackets in 
@). The value of 04 is not yet known. The formula (|l^) clearly displays 
the relation between the NLO and NLOr approximations. One sees that the 
NLOr term ri is summed with the NNLO value of 02. The analogous situation 
occurs for higher order terms. In particular, the 4NL0 term contains, beside 
a4, the contributions from the 3NL0r, NNLOr and NLOr approximations. 

We conclude that within the present accuracy of 0(7q) corrections, the 
perturbative QCD approach fails in the precise determination of the logarith- 
mic slope of the quark jet multiplicity 7i? even at the Z° and can be trusted 
only at much higher energies. 

9 



4 Theory and experiment on multiplicities and 
slopes 

Experimental results on multiplicities and slopes in quark and gluon jets are 
available in the energy range from the T to the Z° 0, |^, |l^ . The CLEO data 
at the T and the OPAL data at the Z^ are fully inclusive (unbiased 
jets). Thus, the jet definitions correspond closely to theory. The DELPHI 
data |T^ at intermediate scales to OPAL and CLEO are based on three jet 
events selected using a jet finder. Parameters of the DELPHI distributions 
are given in terms of the jet hardness scale k, = 2E sm{Q/2), where G is 
the opening angle. For small G, one gets y ~ IhIk/Qq) as mentioned above. 
The data of CLEO and OPAL are presented in terms of the jet energy or 
invariant mass. Leaving aside the question of whether the DELPHI results 
can be directly compared with pQCD and the difference in scales at large 
G, we proceed to an analysis keeping in mind that qualitative trends and 
10-20% accuracy will satisfy us at present. 

The experimental results for r are shown in Fig. 1. The data are seen 
to lie substantially below the predictions. The theoretical results approach 
the experimental values more closely, however, as higher order terms are 
included. The comparatively large value of r2 helps in this regard although 
it also poses some problems for the perturbative treatment of slopes, as 
discussed in section ^. Agreement within 10-20% accuracy is achieved for 
the OPAL value of r at the Z°. Here E ^ AO GeV and k ^ 37 GeV is only 
slightly smaller than E. With 05(37 GeV) ~ 0.14 and 70 ~ 0.517 one obtains 
from d^) r ^ 1.71 for Uf = 3 and r ^ 1.68 for Uf = 4, i.e. much smaller 



values than 2.25 at asymptotics. The computer solution |TT[ leads to perfect 
agreement with the OPAL result r^xp ~ 1.51. The computer solution differs 
from analytical estimates by use of the explicit multiplicities in the integrals, 
with no Taylor series expansion. This fact becomes especially important for 
higher order terms at present values of the expansion parameter 70. The 
limits e~^ and 1 — e~^ of integration and the scale of the running coupling 
constant (however, in the one-loop approximation only) are included. 

According to eq. (|^), the ratio r becomes smaller with decreasing scale 
since 70 increases. The same is observed for the experimental results. Thus 
the qualitative trends agree. The computer solution [|TT| decreases more 



slowly than the experimental data, exceeding the experimental value of r 
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by 15-20% at the T as mentioned in the introduction. This results in a 
noticeable difference between the theoretical (computer) and experimental 
values of r . From Fig. 2 in |TT[ , the slope predicted by the computer solution 
can be estimated to be about 0.096 at the scale of the Z° . From the data 
in Fig. 1, the corresponding result is about 0.174. These values differ from 
the analytic prediction of 0.06 given in section ^ This indicates the strong 
influence of higher order perturbative corrections, even at the Z^, and of 
nonperturbative terms at lower energies which reduce the energy dependence. 
We have also performed a fit of the gluon jet multiplicity. The only 
dependence on energy in (^), (^ comes from as{y) which we constrain to 
agree with the PDG results ||2l|. The energy dependence of the gluon jet 
multiplicity at 3NL0 is given by 



nG) = Ky--'^ exp[2Cy^ + 5G(y)], 



(18) 



with K an overall normalization constant, C = JANc/Po: ^ind 



Sciy) = -p 



2a,C' + %[\n{2y) + 2] 



y 



asC' - ^[H2y) + 1] 
Po 



(19) 

The normalization of the gluon jet multiplicity is arbitrary. It is determined 
by the lower limit of integration yo in 



(ug) oc exp 



i{y')dy' 



yo 



(20) 



which is not defined. We perform a one parameter fit of eq. ( [TSD to the gluon 
jet data in []T^, with K the fitted parameter. Generally good agreement with 
this data is found as shown in Fig. 2. The rate of the multiplicity increase 
is somewhat overestimated by the theory, however, especially by the 3NL0 
(dashed) curve. The power-like terms provided by the nonperturbative region 
(i.e. the truncation of the integration in Vj at e~^ and at 1 — e~^) are not very 
important. They slightly flatten the energy dependence shown in Fig. 2. 

One can write the energy dependence of the quark jet multiplicity anal- 
ogously to (|T8|) as @ 



To 



(21) 



11 



with 

SF{y) = Sciy) + ^ri + —{r, + ^). (22) 

Vy y ^ 

If the normahzation factor K is determined by the fit to the gluon jet data as 
in Fig. 2, the normahzation of the quark jet multiphcity according to formulas 
(H), (H), (|]) is too low by about 2-3 units, as shown in Fig. 2. This is due to 
our failure to describe the ratio r with sufficient accuracy as discussed above. 
More important is the fact that were we to normalize (np) by brute force 
at some energy, we would be unable to describe its slope. This can be seen 
from Fig. 2. The quark jet curves with normalization fixed by the gluon jet 
data are seen to have about the same slope as the data. Thus, if the quark 
jet predictions are multiplied by a constant factor so that the experimental 
and theoretical results for (np) agree at some scale k, the experimental and 
theoretical slopes would differ drastically. In particular, the solid curve at 
NLO (i.e. for 6f = 0) would be much steeper than the experimental data 
and the 3NL0 curve would be steeper yet. Nonperturbative corrections 
would flatten the slopes and yield better agreement between data and theory. 
Nonetheless the slope of (np) is less well determined than that of {nc). It is 
more sensitive to higher order corrections (see section ^. 

We should emphasize that there are ambiguities with the data in Fig. 2. 
These concern the dependence on a jet flnding algorithm and the choice of 
the scale variable k. Qualitatively similar results to those shown in Fig. 2 
are obtained using unbiased gluon and quark jets which do not share these 
ambiguities [^]. For the unbiased jets, the growth of multiplicity with scale 



is more marked than that exhibited by the data in Fig. 2 (see Figs. 1 and 2 
in [^). The 3NL0 prediction (p!8D is found to describe the unbiased gluon 
jet measurements well, including the slope, with a value of as consistent 



with measurement ||21[]. Thus the disagreement between experiment and 
theory seen in Fig. 2 for the gluon jet slope is possibly due to the lack of 
correspondence between the theoretical and experimental deflnitions of the 
jets. Concerning quark jets, it was shown in [^ that eq. (|2T|) describes the 
unbiased quark jet measurements well using a value of Qq about three times 
smaller than that found for unbiased gluon jets. Such a scale is physically 
acceptable and leads to a reasonable value of as ■ In contrast, a flt of 



the quark jet data in Fig. 2 using the technique described in |g2[ yields 
Qo ~ 10~^ GeV, an unphysical value of the scale. 
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The ratio of slopes r'^^^ is less sensitive to higher order corrections and is 
better approximated by the NLOr expression than r, as discussed in section |[ 
According to DELPHI |T^, the experimental value r^^Jp depends very little 

on scale and stays almost constant at about 1.97±0.10, whereas rj^^o ~ 
2.01 - 2.03. Using formulas (|), (|10D, (|ll one gets r^^^LO ~ 1-86 - 1.92 from 



the T to the Z°. Fig. 3 displays the values of r'-^-' in the NLOr and 3NL0r- 
approximations for nj = 4 in comparison to the experimental limits [1^ . 



The measured values of pi^exp ~ ^expf^ range from 0.53 at the T to 0.75 
at the Z°. Its perturbative values according to (|1^) range from 0.85 to 0.9. 



Either next order perturbative or nonperturbative terms are needed |1J] to 
reproduce the big difference di = r^^^ — Vexp ~ 2(1 — pi^exp) observed in 



experiment. The simplest Pade expression (|l^) improves the situation only 
slightly. 

There is a widely held opinion p3|, |2^ that the average multiplicity in 
e~^e~ collisions is well described by pQCD in the NLO-approximation. At the 
same time, one does not usually mention that the NLO-approximation fails 
badly in its prediction for r. However these two facts are tightly connected if 
one analyzes multiplicities in separated quark and gluon jets. As the simplest 
approach, one used to consider the formula rie+e-i^E) = 2nF{E) and, more 
importantly, to fit its energy behavior using the anomalous dimension 7 
whereas ■jp = pi'j should be used. In NLO, their ratio pi equals 1 so that this 
replacement seems justified were it not for the completely wrong (compared 
with experiment) NLO- value of r in the definition of pi = r/r^^\ whereas r^^' 
in the NLO-approximation approximately equals its experimental value r^^X- 
Therefore, one should view this success with some scepticism. This issue is 
also discussed in [^. 

This problem revealed itself in experiment when the DELPHI data on 
the scale dependence of multiplicities in quark and gluon jets became avail- 
able |T^. The quark jet logarithmic slope 7^ was found to be much smaller 
than the gluon slope 7, and their ratio pi was measured to be about 0.50-0.75 
as mentioned above. If one insists on using NLO expressions, it is therefore 
necessary to change (jip) while leaving {np) unchanged. This was done 
in [^ by adding an ad-hoc constant to {rip) which was found to be rather 
large (about 2-3). 

It seems premature to introduce constant offsets to the quark and gluon 
jet multiplicities before other, better motivated possibilities are studied both 
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from the theoretical and experimental sides. Apart from the ambiguities in 
the DELPHI measurements discussed above, there is a rather regular way 
to account for higher order perturbative and nonperturbative corrections 
which leads in the proper direction. This is supported by the computer 



solution [^ of eqs. (1), (2) with quite satisfactory agreement with experiment 
for the value of r at the Z^. The slight difference at low energies results in 
a disagreement on slopes which becomes a stumbling-stone for the whole 
approach and requires very accurate calculations. To predict values of the 
slopes with a precision of several percent one needs computer solutions of 
eqs. (1), (2) with the refinements mentioned at the end of section 0. At 
present, analytical estimates provide accuracy on the order of 10-20%. By 
itself, it is astonishing to obtain any agreement at all given the large value 
of the expansion parameter 70 ~ 0.5 at hand. This should be considered as 
a success of the whole approach. 

The predictions about curvatures and widths, as well as higher moments. 



should also be confronted to experiment (see [^, |2g, |27 



5 Summary 

We have shown that purely perturbative QCD estimates of the anomalous 
dimension 7 and the ratio of average multiplicities in gluon and quark jets 
r are quite reliable at the level of about 10-20%, where the nonperturbative 
corrections become important. The simplest estimates of these corrections 
arising from the limitation of the parton cascade by some virtuality cut-off 
are given. 

For the slope of r (i.e., r ) and for the ratio of the logarithmic slopes 
of average multiplicities pi, the perturbative results become reliable only at 
much larger energies than the Z^. However, the qualitative trends of the 
theoretical corrections demonstrate a tendency to yield better agreement 
with experiment as higher orders are included. In particular, the partial 
compensation of higher order corrections in the ratio of the slopes of average 
multiplicities r^^^ results in a good description of the data by the lowest NLO- 
expression. These conclusions are obtained using corrections of order 7q in 
7 and r as given by the equations for the generating functions. To evaluate 
the slopes at the Z^ and below with better accuracy one should solve the 
equations numerically. 
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In the same approximation we have derived the analytical predictions of 
perturbative QCD for the curvatures of the average multiplicities of quark 
and gluon jets as a function of energy (or of the evolution parameter). Recent 
progress in disentangling quark and gluon jets leads a way to the experimental 
verification of these predictions. 

Appendix 1. The anomalous dimension and 
the ratio of multiplicities 

To obtain the perturbative solution of eqs. (1), (2) we use Taylor series 
expansion (at large y) of multiplicities in the integrals and then differentiate 
both sides of the equations. The arguments of 7o in the integrals are chosen 
to be y + lnx(l — x) in accordance with the dependence of the coupling 
strength on the parton transverse momentum. Finally, we get the following 
equations: 



n, 



ANr 



\\unG — '^G ~ "^unp + 2n^)v4 + 2{unQ — Surip + n^jvs + n^vejj, (23) 



2 r / \ / \ X /// 

rorip = 'JqIug + [unc - n^jvr + [2un(. - n^Jvs - -^n^Vg + 

(wn^ + unp - np)vio - -npVu\, 



(24) 



where u = 2B'~f^ and Vj are the following integrals with the limits of integra- 
tion chosen as e and 1 — e evaluated up to terms 0{e) (see also p9[]): 



V2 



Vi = / Vdx 
ln(l — X 



1 X + X jdx = e, 

2 2 ^ 24 ' 



X 



V3 



-2Vlnx(l -x) 



36 



In^(l-x) 



X 



-2\/(ln2x + ln^(l-a;)) 
413 



dx 



(25) 



67 - Qir'^ 
dx = — \-e\iae, (26) 



^^^^^"108+^^^'''^"^^''^ + ^^' ^^^^ 
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X 



V4 = / [x^ + (1 - xf]dx = - - 2e, 
J o 

13 
I'l — xS^] hixdx = h eil — hie), 



li 



V6 



(28) 
(29) 
(30) 



j [x"^ + {I - xf]\Yi^ xdx = — -£(ln^e-21n£ + 2), 

J^dx = J{l-^)dx=^--^-e, (31) 

(32) 



V7 



V8 



$ In xdx 

15 



h£(l -Ine 



Vg 



/io 
$ln^a;(ix = E{\n^ e - 2\ne + 2) 
8 



Vll 



Vio = 

a; 






7r2 5 1 



ln(l - x)dx = — — - + -e{\ne - 3), 
o 8 2 

111^(1 - x)dx = - - 2C(3) + -^(In^ £ - 2 In £ + 2) 



(33) 
(34) 
(35) 



( means the ^-function, with C(3) = 1.2021. 

Using formulas (|), (|^) and keeping terms up to 0{jq) order on both 
sides of eqs. (p3D, ( ^4]) we get 

ai = vi + ^^ f 1 ^ V4, (36) 



02 



Oi 

Vi 

2 12A/', 



^/ 



!!La _ 1)) („. + 2B) - 1 - -— 

ro / 2 12A^c?^o 



(2ri + 3v5),(37) 



as = a^ia, - Vi + -5) - -55i + (ai + -fi)v2 - -vg - ^{[^2 (^ - 7) - 



^0 



-(oiri + 3i?ri - rf - r2)]v4 + 2 



25 - (2ai + 7B- ri) 



^0 



r-2 = n(v7 -ai 



ri = 2vi — vy 



ri/ 



vsH V6},(3^ 

2\ 



^ 1 

o ,. ) - V2 - Vg Vio 

SNcT-Q ro 



^/ 



- [(ro - 2) (ai + 25) (3v4 - 2) - 6V5] , (40) 



12iV,ro' 
rs = 2a3 - 2(02 - r2)(ai - n) - n[{a^ - nf + ABa^ + 35^] + B[Bi + 2vio 

-3a2 + 5r2] + OsVy + (2ai + 55) vg - -vg + — [(2ai + 35 - ri)vio - -vu] + 55r^.(41) 

2 ro 2 
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The numerical values of Oj, r^ are given in Table |1] for 6 = 0. The correspond- 
ing curves for r using n/ = 4 are shown in Fig. 1. 

The formulas for ai,a2,ri coincide with those in ]§, 0. The expression 



for r2 was obtained in |jTO[ . Use of Taylor series expansion allows energy con- 
servation to be included in the vertices of Feynman graphs. This important 
fact reveals itself in a special form of the integrals V2, vio containing vr^. Just 
such a form of corrections was also obtained in the Lund model [^ when 
taking into account energy conservation by cutting off the fractal triangles 
in the {y, /i;T)-variables. 

Earlier, in the renormalization group approach [^], only part of our ex- 
pression for r2 was obtained (the terms with ri in (|40|) ), which is 6 times 
smaller than the other terms. It is difficult to account for energy conserva- 
tion in that approach. 

The rather large value of r2 compared with ri raised suspicion about the 
nature of the series for r and its convergence, even at the Z^. However, the 
much smaller value of r^ suggests that this is just a numerical fluctuation, 
and that the series is well behaved and can be trusted. This is supported by 



results of the computer calculations |11 



Appendix 2. The slopes 

The slope of the ratio of multiplicities is given by eq. (|^). Using the constants 
in Table |l| for nj = 4 we get 

r' ^ 0.167o'(l + 4.970 + 1.57o + 0(70))- (42) 

The factor in front of the bracket is very small in the available energy region 
with 7o ~ 0.5. The second term is, however, quite large compared to 1 
at present energies, even though it vanishes asymptotically. Thus one can 
hardly trust the perturbative expression for r at the Z^ and can use it for 
qualitative estimates only. 

The perturbative expression for pi is 

pi = 1 - Bri-f^[l + (ai + ri + ^)7o + 
(ZH^ + ain + 3r2 + — + a2 + aj + rl + 51)7^]. (43) 
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The correction to 1 is small and vanishes rapidly as the scale increases. How- 
ever, the subsequent terms in the brackets are not small at present energies 
as seen from their numerical values presented in eq. (|I^) for nj = 4. 

The situation is similar for the ratio of the curvatures (second derivatives) 
of multiplicities, which has the solution 

P2 = l- 2Bn-f^[l + (^ + n + ai - -)7o 

I ^^^^3 , o , 2 , , 2air2 , 2 , R 

r\ ri 

2^'^ 25r,-25a,-^)7o^]. (44) 



ri 2 2 

Numerically (for nj = 4), 

P2 = 1 - 0.1472(1 + 5.2I70 + 2.147^). (45) 

Once again the terms in the brackets are larger or compatible with 1 at 
present values of 70 ~ 0.5. All the above expressions can be trusted quanti- 
tatively only at extremely high energies where 70 < 0.2, i.e. as < 0.02, far 
below the present value ag ~ 0.12. 

Thus we conclude that purely perturbative QCD as described by equa- 
tions (1), (2) for the multiplicities with integration limits equal to and 1 
can hardly be accepted for quantitative estimates of the logarithmic slopes. 
Qualitatively, it predicts a positive slope for r and an excess of the ratio of 
slopes r'-^-' over the ratio r, in accordance with experimental findings. How- 
ever their numerical values are lower than the experimental ones if only terms 
quoted above are taken into account. 
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Figure 1: The ratio of the average muhiphcities in gluon and quark jets 
versus the scale k, in the DLA, NLOr, NNLOr and 3NL0r approximations 
using integration hmits 0-1, and in the 3NL0r approximation using hmits 
e~^ and 1 — e~^ ('^(^) hne), in comparison to data [^, ^ |I^. The theoretical 
results are obtained using Uf = 4. 
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Figure 2: The average multiplicities of gluon (G) and quark (F) jets in the 
NLO and 3NL0 approximations, compared to data [0. For the theory, the 
gluon jet normalization is fit to the data; the normalization of the gluon jet 
curve fixes the normalization of the quark jet curve. The theoretical results 
are obtained using n/ = 4. 
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Figure 3: The ratio of slopes of the average muhiphcities in gluon and quark 
jets, r*^^\ in the NLOr and 3NL0r (with its Fade expression) approximations. 
The theoretical results are obtained using rij- = 4. Upper and lower bounds 



of experimental results |]TB|] are shown by the dash-dotted lines. 
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